Bunching visibility for correlated photons from single GaAs quantum dots 
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We study photon bunching phenomena associated with biexciton-exciton cascade in single GaAs self- 
assembled quantum dots. Experiments carried out with a pulsed excitation source show that significant bunching 
is only detectable at very low excitation, where the typical intensity of photon streams is less than the half of 
their saturation value. Our findings are qualitatively understood with a model which accounts for Poissonian 
statistics in the number of excitons, predicting the height of a bunching peak being determined by the inverse of 
probability of finding more than one exciton. 



I. INTRODUCTION 

Semiconductor quantum dots (QDs) are nonconventional 
light emitters which produce nonclassical photons. Since the 
first demonstration of photon antibunching in single QD pho- 
toluminescence (PL) iH 0], numerous attempts have been 
made to develop efficient single-photon sources based on QDs 
embedded in micro- and nanostructures iH 0, IH |^ 01 • 

The nonclassical nature of photons has commonly been ver- 
ified by measuring coincidence probability for the arrival of 
two photons, utilizing a Humbury, Brown and Twiss setup. In 
the case of single photons, the coincidence histogram shows 
an antibunching dip (negative peak), confirming zero proba- 
bility of finding two photons at the same time. The depth of 
the histogram represents a figure of merit for a single photon 
source. For correlated photons, on the other hand, the his- 
togram would present a bunching peak, suggesting a higher 
probability of finding a particular pair of photons. 

Two photons sequentially emitted with a biexciton-exciton 
cascade are an example of the coiTelated photons iSH fioll . 
and are recognized as a resource of entanglement lUlll . Po- 
larization entangled photons in this framework were demon- 
strated in bulk CuCl fT?]; subsequently the researchers im- 
proved the visibility of coincidence, reaching violation of 
Bell's inequality, by adopting high repetition pulses which 
engender suppression of 'accidental coincidence counts^ ifisl 
UM . Triggered entangled photons have recently been demon- 
strated in semiconductor QDs |[T5[|T6ll . Once again, achieving 
high fidelity has becomed a central issue iUtI [Tsll . and co- 
incidence data are often analyzed after making 'background 
subtraction'. These studies imply that weak excitation yields 
high fidelity, i.e., high bunching peaks with low backgrounds, 
while the physical origin has not been clarified. 

In this work we are studying the excitation power depen- 
dence of photon bunching associated with biexciton-exciton 
cascades. We have found that a significant bunching feature 
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FIG. 1 : Schematic representation of the observation of high bunching 
peaks for correlated biexcitonic (B) and excitonic (X) photons. 



is only present at low excitation, and disappears at high exci- 
tation where the intensity of biexciton/exciton PL approaches 
its saturation level. We demonstrate that a smaller visibility at 
higher excitation is not due to normal high excitation effects, 
such as the onset of emissions from multiple carrier states, a 
wetting layer, or incoherent carrier scatterings. Instead, our 
findings are understood by a purely photon-statistical effect, 
where a bunching peak decreases with increasing probability 
of finding each photon in a pulse, as is illustrated in Fig. [T] 
This is quite surprising because a regulated sequence of cor- 
related two photons does not show any bunching feature, al- 
though such a light source should be favorable for practical 
applications. 



n. EXPERIMENTAL PROCEDURE 

The experiments were performed on GaAs self-assembled 
QDs in an Alj.Gai_j.As {x = 0.26 ± 0.01) barrier grown on 
a (311)A surface by droplet epitaxy ITqI |25i . Atomic force 
microscopy demonstrated the formation of lens-shaped QDs 
of 80 nm in diameter, 10 nm in height, and 5 x 10^ cm^^ in 
density. The QDs were embedded in a two-dimensional pho- 
tonic crystal (PhC) membrane, where air holes normal to the 
surface were regularly made with C^y symmetry. In the center 
of PhC, three missing holes were aiTanged in line, forming L3 
defect cavity. The PhC membranes of 140 nm in thickness. 



2 




n I \ I r 

706 708 710 712 714 716 71! 
Wavelength (nm) 



o 
O 

<u 
o 
c 

0) 
"D 
O 

c 
o 
O 



40- 



20- 



50- 

0- 
100 



1 40 nW 




80 |iW 

U, a 




120 ^w 



20 



40 60 
Time (ns) 



80 



FIG. 2: (color online) A series of microPL spectra for a single GaAs 
QD at various excitation intensities. Photon energies for excitonic 
and biexci tonic lines are presented by • and □, respectively. 



FIG. 3: (color online) Histograms for the coincident arrival of two 
photons with biexcitonic and excitonic transitions with 164 ps time 
bin for three excitation intensities. The integration times were 90, 30, 
and 15 minutes for the top, middle, and bottom panels, respectively. 



204 nm in lattice constant, and 80 nm in the air-hole diame- 
ter, were fabricated by electron beam lithography and reactive 
ion beam etching [21, '22']. Spectral characteristics in the PhC 
membranes are presented elsewhere |i23J. 

For the optical study, we used second-harmonic output of 
an optical parametric oscillator synchronously pumped by a 
mode-locked Ti-sapphire laser. The laser system produced ex- 
citation pulses with 3 ps duration and 76 MHz repetition. The 
wavelength was tuned to be 640 nm, exciting the absorption 
edge of the AlGaAs barrier. 

PL from a single GaAs QD was observed with a confo- 
cal microPL setup, using an objective lens of 4 mm focal 
length and 0.42 in numerical aperture. The PL beam was 
split by a 1:1 beam splitter, with each beam fed in a grating 
spectrometer equipped with a fiber-coupled avalanche photo- 
diode (APD). The spectral window for each APD was around 
0.8 meV, being much smaller than the exciton-biexciton split 
in a PL spectrum (see Fig.|2]i. Then, we set the detection wave- 
length of one APD to the exciton line, and that of another APD 
to the biexciton line. Electric output from APD's was sent to a 
coincidence counter, yielding a start/stop event for the photon 
arrival. PL spectra were also monitored by a charge coupled 
device. All experiments were performed at 8 K. 



ni. EXPERIMENTAL RESULTS AND DISCUSSIONS 

Figure |2]presents the emission spectra of the QD which we 
will be examining in this study. At low excitation, a single 
emission line (•) is present at 7 11. 5 nm (1.742 eV), being as- 
signed to PL from neutral excitons (X). With increasing exci- 
tation intensity, another emission line (□) starts to appear at 
713.1 nm (1.738 eV), at the lower energy side of the X line, 
which is assigned to PL from biexcitons (B). The split from 
the X line to the B line is 3.9 meV, which is typical for GaAs 



QDs grown by droplet epitaxy ll24l|25l|26ll . 

At very high excitation, a broad spectral component 
emerges at 709 nm (1.747 eV) with 5 meV in full width at 
half maximum. The similar behavior is observed in PhC sam- 
ples at high excitation. It has also been found that the en- 
ergy depends systematically on the lattice constant of PhC 
12311 ■ We therefore attribute this emission to the highly-excited 
PL continuum enhanced by the cavity resonance. Note that 
we choose a PhC cavity with a relatively low quality factor 
(Q ^ 300), allowing both B and X lines to be on resonance, 
while retaining their high emission efficiency. For the coinci- 
dence study, we will limit excitation power to below the value 
corresponding to the onset of the cavity emission. 

The autocorrelation trace of the X line (not shown) presents 
an antibunching dip at coinciding times 39.9 ns in Fig. [3]), 
demonstrating its single-photon behavior. Coincidence his- 
tograms for cross-correlations between the B and X photons 
are presented in Fig. [3] At low excitation with 40 /iW, the his- 
togram shows a high central peak following sequential back- 
grounds. The presence of a high coincidence peak confirms 
that two photons were generated with a single radiative cas- 
cade. The sequential background is due to coincidence counts 
between photons in temporally separated pulses, which were 
emitted synchronously with the laser source. The bunching 
visibility, i.e., the ratio of a central coincidence peak to back- 
ground side peaks, is evaluated to be 2.7 (±0.1). 

When excitation intensity increases to 80 /iW, we find a 
remarkable reduction in the relative height of a bunching 
peak, while the histogram shows a higher signal-to-noise ra- 
tio, reflecting greater counting rates. The visibility in this 
case is evaluated to be 1.3 (±0.05). For excitation intensity 
at 120 fiW, the coincidence peak further decreases to 1.06 
(±0.02), as shown in the bottom panel of Fig.|3] We should 
note that negligible bunching peak was also observed in our 
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prev ious work, which studied QDs without PhC processing 
02711 . suggesting that this is a general feature related to the 
high excitation regime. 

We are analyzing the power dependence of coincidence his- 
tograms in terms of photon number statistics: Preliminary the- 
oretical work was reported in ll28ll . In this model, we assume 
that the probability of finding B and X photons is simply de- 
termined by the number of photoinjected excitons, obeying 
the following Poissonian distribution function. 
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where N is the mean number of excitons. Since an X photon 
is generated when the number of excitons is more than one, 
and both B and X photons are generated when more than two 
excitons are present initially, the probability of finding an X 
(B) photon in a pulse, Px (Pb), is given by; 

Px = Y.^Nin)^l- exp(-iV), (2) 

n>l 

^B = II^iv(") = l-exp(-7V)(l + 7V), (3) 

n>2 

where we used a relation, J2n>o ^i^) ^ ^■ 

A coincidence peak, g^^^^ (0), is given by a joint probability 
for both B and X photons being present in a pulse, and of 
being counted by two detectors, whose counting yield is ijb 
and rjx, respectively. Thus we find. 
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n>2 



P{n)riBrix = PbVbiIx- 



(4) 



Note that the parameter of rjB and rjx accounts for every effect 
which causes a counting loss, including a finite efficiency of 
photon extraction and photon detection. 

(2) 

A coincidence background, g^Q, is written by a product of 
counting probability for B and X photons, each belonging to 
an uncorrelated pulse. Thus, we find. 
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Figure |4] illustrates expected coincidence histograms for 
various values of N. It shows a central bunching peak follow- 
ing constant side peaks. The relative height of the central peak 
is found to depend on N, and it is larger for smaller numbers 
of N. The bunching visibility is, therefore, expressed by. 
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(6) 



The above equation suggests that a bunching peak becomes 
higher as the probability of finding a dark pulse emitted from 
"zero" excitons increases, as is schematically shown in Fig.[T] 
Note that the above simple expression is also obtained by rig- 
orous formulation based on the quantum regression theorem, 
as is developed in Appendix. 

The bunching visibility sharply decreases as the mean num- 
ber of excitons increases, as shown by the solid line in 
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FIG. 4: Coincidence probability for correlated photon pairs for vari- 



ous values of mean exciton number, A'^, for 773 = rjx = 10 
horizontal (temporal) axis is scaled by pump repetition, Td- 
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Fig. HJa). The same dependence of probability of finding B 
and X photons, which is proportional to the relevant PL in- 
tensity, is presented in Fig. lUb). Through the comparison 
between Figs. [Sja) and |3b), we find that a very low injec- 
tion of excitons is necessary to realize a significant bunching 
peak. As a guideline, we plot a vertical line which intersects 
the visibility curve at g'-^^ (0) /g^Q = 2, the value being avail- 
able with classical incoherent light. It is shown that to achieve 
high bunching visibility beyond classical criteria, we should 
keep the number of excitons lower than ln2 = 0.69. Note 
that this number is much smaller than two, which is a rough 
index for a biexciton to be present in a QD. Correspondingly, 
the intensity of the B (X) line should be less than 0.153 (0.5) 
of its saturation intensity. 

Let us compare the above theoretical prediction with exper- 
imental data. For simplicity, we assume the number of exci- 
tons being proportional to the excitation intensity. Then, we 
evaluate the efficiency of photoinjection through the fit of the 
X and B intensities to the power dependence of Eqs. |2] and 
[3] respectively. The result of this fit is plotted by diamonds 
(X) and circles (B) in Fig. lUb). These PL intensities agree 
with this model. With the use of this injection efficiency, we 
can estimate the number of excitons for each coincidence his- 
togram. The observed data for the coincidence visibility as 
a function of exciton number are finally plotted by solid cir- 
cles in Fig. |3a). We find that the experimental trend, i.e., a 
steep reduction in bunching visibility is reproduced quite well 
by this model, while the observed visibility is significantly 
smaller than the theoretical ones. 

One likely reason for this discrepancy is the influence of 
emission from multiple carrier states (or a wetting layer). The 
onset of these incoherent signals results in the emission of un- 
correlated photons, thus the reduction of bunching visibility. 
Moreover, in the above treatment, we assumed that excitons 
followed Poissonian statistics, and their mean number was 
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5: (color online) (a) Relative height of bunching peaks, 



J2) 



as a function of mean numbers of excitons. The solid 



circles show the bunching visibilities that are obtained experimen- 
tally, and the solid line shows the theoretical dependence which fol- 
lows Eq. |6] . For comparison, the same power dependence of the 
intensity of exciton and biexciton emissions is plotted in (b). 



proportional to the excitation power However, neither hy- 
pothesis is, in fact, relevant at high excitation when the exciton 
number approaches the maximum number of carriers which 
can be occupied by a QD. At such high excitation both B and 
X photons are likely more saturated than in this model, thus 
the actual visibility becomes lower than that of Fig. 13 a). 



IV. CONCLUSIONS 

In conclusion, we have demonstrated the fidelity of bunch- 
ing statistics associated with biexciton-exciton cascades de- 
pends on excitation intensity, and significant bunching feature 
only appears at very low excitation such that iV <C 2. This 
situation is in stark contrast to the experimental condition for 
characterizing a single-photon emitter, where an antibunching 
dip is normally measured at sufficiently high excitation, real- 
izing regulated single-photon pulses, and high counting rates. 
On the other hand, the usage of dilute photon pulses is es- 
sential for characterizing the bunching feature for correlated 
photons, although it is often difficult to take coincidence mea- 
surements at sufficiently low excitation. 

In atomic physics a model for a three-level system un- 
der low continuous-wave excitation describing the power de- 
pendence of photon bunching has been derived |29]. Here, 
we have analyzed photon statistics between biexcitonic and 
excitonic recombinations with pulsed excitation. Note that 



the similar power sensitive behavior should be involved in 
polarization-resolved cross-correlation experiments. In this 
case the reduction of coincidence fidelity is induced by sev- 
eral microscopic effects such as a fine level split of excitons, 
spin relaxation, and the onset of unwanted emissions, as well 
as the photon saturation effect which we have studied. 
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APPENDIX: SECOND-ORDER CORRELATION FUNCTION 
FOR SEQUENTIAL PULSED EXCITATION 

We will be deriving the expression of the interbeam second- 
order correlation function associated with biexciton-exciton 
cascades. Let us assume that QDs are excited by short optical 
pulses, and the pulse interval, T^, is sufficiently longer than 
the relaxation times of QDs. We restrict ourselves to the three 
levels consisting of the biexciton |2), exciton and vacuum 
|0) states, as depicted in the inset of Fig. |6l The presence of 
higher excited states influences the initial population of |2), 
followed by the fast energy relaxation after initial photoinjec- 
tion. The population dynamics is characterized by a set of rate 
equations, 

dp22{t)/dt=-A2P22it), (A.l) 

dpii{t)/dt = A2P22it) - Aipiiit). (A.2) 

where pu denotes the diagonal density matrix element of the 
i exciton level, and A2 (Ai) is a biexciton-exciton (exciton- 
vacuum) transition rate. The above equations have the general 
solutions. 



Piiit) = - 



P22{t) = P22i0)e-^'\ 
^2^22(0) -A.t , r /n\ ^ 



(A.3) 
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The second-order correlation function is presented in terms 
of dipole projection operators, tt2 — |1)(2|, vri = |1)(0|, and 
their conjugates |^29l]. 



gg\t,t + T) = 



{TTl{t)TTl{t + T)TT,{t + T)TT2{t)) 
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(A.5) 



The quantum regression theorem allows to express double- 
time expectation values in the right hand side of Eq. IA.5I in 
terms of single expectation values. Thus, we obtain. 
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FIG. 6: (color online) (a) Calculated second-order correlation func- 
tion for Ai — A2/2 — Td/10, and Poissonian distributed exciton 
population with N = 1, which corresponds to /922(0) — 0.264 and 
pii(O) ~ 0.6321, according to Eqs. 2 and 3. (b) Measured coinci- 
dence histogram at 80 /iW. 



In the experiments, the coincidence counts are acquired 
over a long integration time. Tint ^d)- The time-integrated 
coincidence counts are, therefore, given by. 
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h 
h 



A2 



(A.7) 
(A.8) 
(A.9) 

(A. 10) 
(A. 11) 



Note that the value of Ii (I2) is proportional to the averaged 
intensity of the biexciton-exciton (exciton-vacuum) transi- 



tion. Equation |A.9| presents that G^"^^ (t) = for r < 0, and 
it decays with the exciton decay rate for r > 0. In addition, 
the shape of G^^-* (t) is independent of the initial population. 

Next, we consider photon coincidence events between un- 
corrected pulses, which cause background side peaks in co- 
incidence histograms. In this case, double-time expectation 
values in Eq. IA.5I are decomposed in the product of two ex- 
pectation values; 

(4 {t)nl {t + Td + r)^i [t + Td + t)tt2 (t)) 

- {4it)Mt)){4it + r)Mt + r)). (A.12) 
Thus, we obtain the time-integrated coincidence counts given 

by. 
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(A. 13) 



where 



{A2 ~ ^l)/9ll(0)p22(0) + ^2/922(0)2 



P22{0? 

2{A2-A,) 



(A. 14) 
(A.15) 



Equation IA.13I shows that each side peak rises with the 
biexciton decay rate, and it decays dominantly with the exci- 
ton decay rate. Calculation results for the second-order coiTe- 
lation functions are presented in Fig. 6(a), together with an ex- 
perimental coincidence curve in Fig. 6(b). Asymmetric shape 
in both central and side peak is well reproduced by the sim- 
ulation. However, we can only deal with the area of G'^^ (r) 
for quantitative analysis, because of the fast recombination of 
GaAs QDs, whose time scale is similar with that of the in- 
strumental response function. In this case, bunching visibility 
defined by Eq. |6]becomes. 



CooG^Bcir + Wr Pii(0)+P22(0) 



(A. 16) 



Note that Eq. IA.16I is equivalent to Eq. |6l because the value 
of {/5ii(0) + P22(0)} is nothing but the probability of finding 
more than one exciton in a QD. 
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